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ABSTRACT 

Extended material at large radii surrounding a supernova can result in a double-peaked light curve. This 
occurs when the material is sufficiently massive that the supernova shock continues to propagate into it and 
sufficiently extended that it produces a bright first peak. Such material can be the leftover, low-mass envelope 
of a star that has been highly stripped, the mass associated with a wind, or perhaps mass surrounding the 
progenitor due to some type of pre-explosion activity. I summarize the conditions necessary for such a light 
curve to occur, describe what can be learned about the extended material from the light curve shape, and 
provide an analytic model for fitting the first peak in these double-peaked supernovae. This is applied to the 
specific case of a Type Ic super-luminous supernova, LSQ14bdq. The mass in the extended material around 
this explosion’s progenitor is measured to be small, ~ 0.3 — 0.5 Mq. The radius of this material can be 
~ 500 — 5,000 but it is difficult to constrain due to a degeneracy between radius and the supernova’s 
energy. In the future, spectra taken during the first peak will be important for measuring the velocity and 
composition of the extended material so that this degeneracy can be overcome. 

Subject headings: supernovae: general — supernovae: individual (SN 1993J, LSQ14bdq) 


1 . INTRODUCTION 

Observations of supernovae (SNe) during the first few days 
of their rising light curves provide valuable information about 
their progenitors and the circumstellar environment of the ex¬ 
plosion (Piro&Nakar 2013). Although historically it has 
been difficult to catch SNe at such early moments, current 
and forthcoming wide-field surveys make this an ideal time to 
put more focus on early light curves. 

One of the exciting results from such observations is the 
discovery of double-peaked SNe, where the first peak lasts 
about a day, the second peak lasts a couple weeks, and the 
luminosity in each peak is relatively similar. It is important 
to note that single bands at short wavelengths can also show 
double peaks (for example, //-band for SN 1987A) due to the 
temperature evolution of the ejecta (first heated by a shock and 
then later heated by radioactive nickel), and this is generally 
expected for a normal massive star (Nakar & Sari 2010). This 
type of double-peaked light curve is distinctly different than 
the cases I focus on here. The key difference is that when two 
peaks are seen in R- or /-band it points to more unique stellar 
structure (Nakar & Piro 2014). 

Perhaps most famous among the double-peaked light 
curves have been the subclass of SNe lib that exhibit this fea¬ 
ture, such as SN 1993J, 201 ldh, and 2013df (Wheeler et al. 
1993; Arcavi et al. 2011; VanDyketal. 2014). It is now 
generally accepted that their first peak comes from the 
presence of low-mass (~ 0.01 — 0.1 Mq), extended (~ 
10 13 cm) material (Woosley et al. 1994; Bersten et al. 2012; 
Nakar & Piro 2014). This unique structure is also reflected 
in pre-explosion imaging, which identified the progenitors as 
yellow supergiants (Aldering et al. 1994; Maund et al. 2011; 
Van Dyk et al. 2011,2014). 

In principle, extended material can be present around a 
wider variety of different SN progenitors, and thus it is use¬ 
ful to develop a flexible theory that can be utilized to study 
all of these cases. Most recently, LSQ14bdq was discovered 
to be a double-peaked super-luminous SN Ic where the first 
peak lasted ~ 15 days (Nicholl et al. 2015). Although such 


an energetic explosion is clearly in a different regime than the 
previously mentioned SNe lib, similar physics should apply. 
Previous work on double-peaked SNe has focused on detailed 
numerical models (Woosley et al. 1994; Bersten et al. 2012), 
which are difficult to use for investigating a wide parameter 
space, or analytic scalings (Nakar & Piro 2014), which do 
not provide time-dependent light curves for direct compari¬ 
son with observations. In the present work I bridge this gap 
by developing a time-dependent light curve models that can 
be fit to the observations. 

In Section 2, I summarize the condition necessary for a 
double-peaked SN light curve and then present simple scal¬ 
ings connecting the first peak light curve to the properties of 
extended material, which have been utilized previously in the 
literature. These scaling are used in Section 3 to develop an 
analytic, time-dependent model of the first peak. In Section 
4, 1 utilize this model to constrain the properties of extended 
material in LSQ14bdq. In Section 5, 1 summarize the conclu¬ 
sions and discuss future work. 

2 . CONDITIONS NEEDED FOR A DOUBLE PEAK AND 
ANALYTIC SCALINGS FOR THE FIRST PEAK 

I begin by describing the conditions needed to produce a 
double-peaked SNe and summarizing the main features of 
their light curves. This relies heavily on the analytic results 
of Nakar & Piro (2014), but I provide them here since they 
are used extensively in the next section. The basic set up is an 
explosion of a core with mass M c and radius R c , which is sur¬ 
rounded by extended material with mass M e and radius II,.. 
As the SN shock powered with an energy E sn passes from the 
core to the extended material, it transitions to a new velocity 
v e . This continued propagation of the shock can only occur if 
M e is sufficiently optically thick such that 

AjrR^ q 

M e > -^—= 6 x v^R^Mq, (1) 

K V e 

where the opacity is K 0.34 = k/ 0.34 cm 2 g -1 , as appro¬ 
priate for electron scattering of solar composition material. 
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vg = u e /10 9 cms -1 , and i?i 3 = f? e /10 13 cm. The veloc¬ 
ity in the extended material estimated considering the reverse 
shock that occurs at the abrupt density drop and energy con¬ 
servation, resulting in 


density U = 3E/4:nR(t) 3 by 

1 4 /1\ 1 4 t; 

TdS = -dU + -Ud - = -dU + —Udt . (10) 

P 3 \pj p pR(t ) 


v e 


2 x 10 9 El i s 


\m q J 


M e 

0.01 Mq 


-0.15 

cms -1 , 

( 2 ) 


where E 5 1 = E sn / 10 51 erg. This implies that only an amount 
of energy 


E e 


4 x 10 49 E 51 


M c \ ~ 0 - 7 ( M e 


0.7 

erg, (3) 


is passed into the extended material. 

The shock-heated extended material subsequently expands 
and cools, producing the first peak. The timescale for this 
emission to reach its maximum luminosity is set by equating 
the expansion and diffusion timescales (Arnett 1982), which 
is roughly estimated as 


t 


p 


f M 6 k\ 1/2 
V v eC ) 


and the peak luminosity is 


( 4 ) 


The entropy decreases due to radiative cooling (I assume no 
radioactive heating in the extended material). Heat balance is 
then given by 


T—M e = -L 
dt 


( 11 ) 


resulting in 


47rf?(f) 3 

3 


dU 

dt 


m 


4 U 


= -L 


( 12 ) 


This expression can be written in a simpler form by substitut¬ 
ing E for U and performing some algebra, 


dE E 
dt t + t e 


( 13 ) 


where t e = R e /v e is the expansion timescale. 
The radiative luminosity is approximated as 


L 


(t + t e )E 

T 1 

z p 


(14) 


L 


E e R e 

V e tl 


( 5 ) 


Equation (4) has the correct scalings, but it is missing a pref¬ 
actor that is necessary for quantitative comparisons with ob¬ 
servations. For this reason, Nakar & Piro (2014) fit to numer¬ 
ical work (Woosley et al. 1994; Bersten et al. 2012), resulting 
in 


tp ~ 0-9 k 0^.34,^51 ^ 


0.17 


Me 


M r . 


m q j v°-° 1M © 


0.57 


day, (6) 


and 


Lp KS 2 X 10 43 K 0 34^51 f?!3 ( —— 


Me 


0.01 Mq 


Mq 

-0.3 


-0.7 


ergs 


( 7 ) 


Although Equations (6) and (7) can be utilized to constrain 
the values of M e and R e from observations, it would also be 
useful to have time-dependent light curves for more detailed 
comparisons as well. This motivates the model presented in 
the next section. 


3. ONE-ZONE MODEL FOR THE FIRST PEAK 

A time-dependent light curve of the first peak is derived 
using a one-zone model where the entirety of the extended 
material is described with characteristic values. The radius of 
the extended material expands as 

R(t) = R e + Vet, (8) 

so that it has a characteristic density given by 

pit) = 3Me/^R{t) 3 . (9) 

Following the first law of thermodynamics, the variation of 
specific entropy is related to the characteristic internal energy 


as is expected from numerous previous analytic works (e.g., 
Arnett 1982; Popov 1993). With this choice of the radiative 
luminosity. Equation (13) is solved analytically, resulting in 


L{t) 


i 2 

Lp 


exp 


tit + 2 t e ) 
2 t 2 p 


( 15 ) 


This correctly gives L « t e E e /tp at early times, consistent 
with Equation (5). The mains limitations of this approach 
is that it does not include corrections to the luminosity and 
temperature during the rise that are dependent on the den¬ 
sity structure of the stellar profile. For a red supergiant struc¬ 
ture, L oc f -0 ' 17 and T e g oc i -056 (for times later than the 
shock crossing time, Nakar & Sari 2010), while, assuming a 
black body spectrum, the model here gives L ~ constant 
and T e ff oc A °' 5 . Since this mostly impacts the earliest parts 
of the rise where data is relatively sparse, ignoring these de¬ 
tails is a good compromise because the main point of the 
model is to capture the peak and fall afterwards. If necessary 
though, these scalings can be used to improve the solutions 
for t < t p , and this would be especially important for bluer 
bands where the temperature evolution makes a larger impact. 
At late times, the luminosity falls as ~ exp(— t 2 /2t p ). This 
solution is not accurate if recombination in the extended mate¬ 
rial begins. In all the cases used in this paper, the light curves 
begin to rise due to radioactive nickel heating before suffi¬ 
ciently low temperatures are reached, and thus recombination 
can be ignored. 

I provide examples of the range of light curves using Equa¬ 
tion (15) in Figure 1. The E-band absolute magnitude is cal¬ 
culated assuming black body emission. A cut off in the light 
curve is added for times less R e /c to replicate light travel ef¬ 
fects. Data from SN 1993J are also shown as an example 
of a prototypical first peak from an SN lib (Wheeler et al. 
1993). The green curve best fits this data, and uses values 
of Me = 0.03 Mq, Me = 2.5 Mq, R e = 4 x 10 13 cm and 
E sn = 0.6 x 10 31 erg. This core mass was chosen to be 
similar to the values typically found in Woosley et al. (1994) 
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Figure 1. Example first peak light curves found from Equation (15). Black 
data points are for SN 1993J (Barbon et al. 1995; Nicholl et al. 2015). The 
green solid curve is the best fit model, which uses M e = 0.03 Mg, M c = 
2.5 Mq, R e = 4 X 10 13 cm and E an = 0.6 X 10 51 erg. The opacity 
is electron scattering for solar-composition material. The other curves use 
all the same parameters, except in each case one parameter is changed. For 
the red dashed curve M e = 0.1 Mq, for the blue dotted curve R e = 5 X 
10 12 cm. and for the purple dot-dashed curve E sn = 2.5 X 10 51 erg. See 
the text for a further discussion of the impact each of these variables have on 
their respective light curves. 

and Bersten et al. (2012). The other curves show examples 
of when one of these values is changed to provide some idea 
about how strongly dependent the light curve is on these pa¬ 
rameters. In particular, M e and E sn impact the width of the 
light curve, as can be seen in Equation (6) where these param¬ 
eters have the strongest dependencies. On the other hand, R e 
and E sn impact the peak luminosity. Typically R e can vary 
for a given stellar progenitors much more (just think about 
the difference between a Wolf-Rayet star and a red super¬ 
giant) than E sn (which is usually around ~ 10 51 erg). This 
is the reason light curves from shock cooling are so useful for 
constraining radii (as emphasized in Piro & Nakar 2013). 

Another aspect to note is that the best fit extended mass of 
M e = 0.03 Mq is much less than values of ~ 0.1 Mq that is 
typically presented by Woosley et al. (1994) and Bersten et al. 
(2012) for SNe lib. This is because in fact the first peak is 
most sensitive to only the mass near the maximum radius and 
not the total amount of hydrogen present (see the more de¬ 
tailed discussion in Nakar & Piro 2014, and in particular their 
Fig. 2, which explicitly shows how the mass measured by the 
first peak compares with the total hydrogen shell mass). So 
while the total hydrogen can indeed be ~ 0.1 Mq to produce 
a realistic, hydrostatic model, the first peak itself only can be 
utilized to measure M e = 0.03 Mq. With this model in hand, 
I now turn to a very different kind of double-peaked SN to 
explore if this same model can be applied. 

4. APPLICATION TO LSQ14BDQ 

Super-luminous SNe Ic are hydrogen/helium-poor explo¬ 
sions that reach peak absolute magnitudes of brighter than 


—21 (Quimby et al. 2011; Gal-Yam 2012). Although they are 
relatively rare (less than ~ 1% of SNe, Quimby et al. 2013; 
McCrum et al. 2015), their huge luminosities make them ob¬ 
servable at cosmological distances. Nevertheless, their un¬ 
derlying power source is still unclear. Unlike super-luminous 
SNe II, they do not show a clear signature of interaction 
with circumstellar material. So although this in principle is 
still a potentially viable model for generating their luminos¬ 
ity, other scenarios like spin down of a newly born magnetar 
(Kasen & Bildsten 2010; Woosley 2010) and fallback accre¬ 
tion onto a black hole within the SN (Dexter & Kasen 2013) 
are also under consideration. There is considerable diversity 
though, and some super-luminous SNe Ic are even candidates 
for pair-instability SNe (Gal-Yam et al. 2009; Young et al. 
2010) where a massive core (greater than ~ 65 Mq) under¬ 
goes complete thermonuclear disruption (Heger & Woosley 
2002 ). 

Most recently, LSQ14bdq is a super-luminous SN Ic 
that was discovered to have a double-peaked light curve 
(Nicholl et al. 2015). Previously, SN 2006oz was another 
super-luminous SN with an early peak (or at least an early 
plateau, Leloudas et al. 2012), but the early detection of the 
rise of LSQ14bdq makes it an especially useful case for ap¬ 
plying this theory. Although there is some uncertainty in the 
exact time of explosion, the first peak reaches a maximum 
luminosity at ~ 7 days after the explosion and a minimum 
occurred ~ 15 days after explosion. The brighter and broader 
main peak occurs ~ 50 days after explosion. Nicholl et al. 
(2015) disfavor having the first peak powered by radioactive 
nickel. The problem is that its luminosity would require more 
nickel than the total mass needed to explain the ~ 15 day 
timescale width (a typical SN la lasts just a little longer but is 
considerably dimmer, showing how a radioactively-powered 
model is difficult to reconcile with the first peak). 

Nicholl et al. (2015) instead favor cooling emission of 
the shock-heated surface layers of the exploding star, sim¬ 
ilar to the models studied by Nakar & Sari (2010) and 
Rabinak & Waxman (2011). From this they infer that the ra¬ 
dius of the exploding star must be large, ~ 500 Rq. While 
such emission can indeed produce a double-peaked light 
curve, as was observed for SN 1987A, this was in /J-band. 
LSQ14bdq was instead observed in U-band at early times. 
In redder bands it is difficult to make a double-peaked light 
curve with the normal structure of a red or blue supergiant be¬ 
cause as the ejecta expands and cools, it gets brighter in red¬ 
der bands even if it becoming less luminous bolometrically. 
Therefore if the first peak is simply shock cooling it would not 
be as tightly peaked. So even if the radius that Nicholl et al. 
(2015) infer may be roughly correct, the amount of mass be¬ 
ing shocked must be much lower than what they consider. 
Instead it should be like the extended material I have been 
discussing here. This event therefore represents an exciting 
new opportunity to study extended material around a totally 
different type of SN. 

Unlike typical SNe, for which E sn ~ 10 51 erg is a resum- 
able assumption, it is much less clear what to use for the 
case of a super luminous-SN. This problem is compounded 
because there are degeneracies between E sn , M e , and R e , 
and only (/-band is available for LSQ14bdq for the first peak. 
For these reasons, I consider both the best fits to the data 
and then try constraining certain parameters to make addi¬ 
tional fits. This helps highlight which values are most cer¬ 
tain and what range is reasonable. The first comparison is 
in Figure 2, where I look at the extended model versus the 
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Figure 2. Example fits to the first peak of LSQ14bdq using three differ¬ 
ent models. The blue and red curves are for stellar models with radii of 
100 Rq and 500 Rq , respectively, similar to those presented in Nicholl et al. 
(2015) using the work of Rabinak & Waxman (2011). These roughly match 
the rise, but are not able to reproduce the drop that begins after ~ 7 days. 
The green curve shows an extended material model with M e = 0.30 Mq. 
M c = 30 Mq, R e = 5,000 Rq, and E a n = 8.8 x 10 51 erg. 

shock cooling of a blue or red supergiant. These use the ana¬ 
lytic scalings from Rabinak & Waxman (2011) with values of 
fZ* = 100 Rq, M* = 60.7 Mg, E sn = 1.52 x 10 53 erg, and 
f p = 0.05 for the blue supergiant model and fZ* = 500 Rq, 
M» = 29.7 Mq, E sn = 1.78 x 10 52 erg, and f p = 0.05 
for the red supergiant model (f p is a density parameter that 
only weakly impacts the light curves, see Rabinak & Waxman 
2011 for a more detailed discussion). These specific val¬ 
ues are chosen, and Rabinak & Waxman (201 1) is used rather 
than Nakar & Sari (2010), to provide results similar to those 
presented in Nicholl et al. (2015). The main conclusion to no¬ 
tice is that these standard stellar models cannot provide a suf¬ 
ficiently quick drop in comparison to the observed first peak. 

In contrast, I also plot in Figure 2 an extended material 
model with M e = 0.30 Mq, M c = 30 Mq, R e = 5,000 Rq, 
and E sn = 8.8 x 10 51 erg. The opacity is assumed to be 
k = 0.2 cm 2 g _1 , since this event showed no signs of hydro¬ 
gen, but the lack of spectra during the first peak mean that 
the exact opacity is uncertain. This core mass is chosen to 
be similar to the value Nicholl et al. (2015) find by fitting the 
second, bright peak of LSQ14bdq with a magnetar model, and 
this value is used throughout the rest of this paper for the ex¬ 
tended material models. Since M c just sets the conversion 
between E sn and E e via Equation (3), using a different value 
of M' c simples implies a new energy of 

E' sn = (M' c /M c )°- 7 E sn . (16) 

Thus, for this model if I instead used M c = 10 Mq, the corre¬ 
sponding energy would need be E sn = 4.1 x 10 53 erg to pro¬ 
vide a similar fit to the data. Figure 2 shows that this extended 
material model provides a much better description of the first 
peak, but the physical picture is now qualitatively very differ- 
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Figure 3. A diverse set of extended material models that can adequately 
explain the first peak of LSQ14bdq. The specific values used in these models 
are summarized in Table 1. 


Table 1 

Extended Material Models for LSQ14bdq 


Me 

(AT©) 

Re 

(Rq) 

Esn 

(10 51 erg) 

E e 

(10 51 erg) 

V e 

(kms -1 ) 

0.51 

500 

54.9 

3.4 

28,000 

0.39 

1,000 

35.0 

1.8 

23.000 

0.30 

3,000 

15.2 

0.66 

16,000 

0.30 

5.000 

8.8 

0.39 

12,000 


All models use M c = 30 Mq and k = 0.2 cm 2 g 1 . Other values for M c 
change the associated E su , which can be inferred using Equation (16). 


ent than what is usually considered. Instead of just a simple 
massive star, there is the massive star surrounded by extended, 
low mass material. 

Unfortunately, the extended model I use in Figure 2 does 
not uniquely fit the observations. As can be seen from the 
strong dependence on both E sn and R e in Equation (7), there 
are degeneracies between these two parameters. For more 
typical core-collapse SNe, assuming E sn ~ 10 51 erg would 
be reasonable, but this cannot be done for a super-luminous 
SN. To highlight how much uncertainty this introduces, in 
Figure 3 I plot a number of different solutions that can fit 
the data, the associated parameters of which are summarized 
in Table 1. This shows that the best constrained parame¬ 
ter is M e , which is set by the width of the first peak. This 
should be expected because the strongest dependence for t p 
in Equation (6) is on M e . This robustly shows that M e is 
relatively low, in the range of ~ 0.3 — 0.5 Mq. One note 
of caution is that this only represents the mass that is mostly 
at R e . Similar to models of SNe lib (Woosley et al. 1994; 
Bersten et al. 2012), if this extended material is the configura¬ 
tion of a hydrostatic envelope, the total mass in the envelope 
may be greater than what is measured with M e . 

On the other hand, there are a wide range of values for E sn 
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and R e that are possible. In particular, R e needs to be large 
andean vary from ~ 500—5,000 Rq. The lower value is simi¬ 
lar to Nicholl et al. (2015), but the larger values are somewhat 
greater than what one might consider for a standard super¬ 
giant model. Above R e ~ 5,000 Rq is not possible because 
the envelope becomes too diffuse and the shock cannot con¬ 
tinue propagating as given by Equation (1). A key way to 
discriminate between these possible solutions would be spec¬ 
tra during the first peak to estimate v e . A larger v, would 
signify a smaller R e , and vice versa. I include associated 
values of v e in Table 1 for reference. Such spectra would 
also measure the composition of the extended material, which 
is also unconstrained for the models I calculate (here I sim¬ 
ply use k = 0.2 cm 2 g -1 , motivated by the lack of hydrogen 
and helium seen at later times). One might think that the fact 
LSQ14bdq is super-luminous, then E sn is relatively large and 
(consulting Table 1) R e relatively small. Unfortunately, this 
does not have to be the case because the high luminosity of the 
second, main peak could be powered by delayed energy injec¬ 
tion (such as in the magnetar model, Kasen & Bildsten 2010; 
Woosley 2010) and thus the energy associated with the first 
peak could be much lower like a typical SN. Again, spectra 
and velocity measurements during the first peak would help 
remove this uncertainty. Mapping the time-dependent veloc¬ 
ity evolution better might also help probe the interior structure 
of the star (Piro & Morozova 2014). 

One deficiency in utilizing these theoretical light curves 
is that they do not constrain the distribution of mass within 
the extended material. This makes it difficult to discern be¬ 
tween an extended shell that is left over from mass loss (which 
would have p ~ r~ 3 or p ~ r _3//2 depending on whether the 
material was radiative or convective, respectively, although a 
large radius would favor a convective profile), a stellar wind 
(which would have p ~ r~ 2 ), or even some other non-trivial 
mass distribution due to pre-explosion activity of the progen¬ 
itor. More detailed numerical investigations of the first peak 
light curve would provide more insight into how well these 
properties can be constrained. 

5. CONCLUSIONS AND DISCUSSION 

The discovery of double-peaked SN light curves is an excit¬ 
ing development that promises to provide new and unique in¬ 
formation about the progenitors of SNe (Nakar & Piro 2014). 
LSQ14bdq has provided an opportunity where the theory for 
such light curves can be applied to a super-luminous SN Ic, 
which can help shed light on the many uncertainties that sur¬ 
round these relatively rare events. Motivated by this, I de¬ 
velop an analytic model to constrain the mass and radius as¬ 
sociated with the first peak in this and other double-peaked 
SNe. This is summarized by Equation (15), which can used 
along with Equations (3) and (6) for setting the initial en¬ 
ergy and diffusion time of the extended material, respectively. 
The main conclusion for LSQ14bdq is that the extended ma¬ 
terial is relatively low mass (~ 0.3 — 0.5 Mq) in compari¬ 
son to the bulk mass associated with the SN. The radius is 
more difficult to constrain, although it must be relatively large 
(~ 500 — 5,000 Rq) to explain the high luminosity of the first 
peak. 

In the future, continued discoveries and studies of double- 
peaked SNe will help inform us about the diversity of progen¬ 
itors that can give rise to SNe. In most cases, the first peak 
will likely not be as broad as for LSQ14bdq, and thus will 


require cadences of less than a day to be properly resolved 
(as has been seen for the double-peaked SNe lib). Short ca¬ 
dence surveys like Zwicky Transient Facility (ZTF, Law et al. 
2009) will be important for such studies. A wide-held UV 
survey (Sagiv et al. 2014) would be especially well-suited for 
this work, since the temperatures in this first peak are rela¬ 
tively high, ~ 10,000 — 20,000 K. The key will be catch¬ 
ing the first peak sufficiently early, so that the time of explo¬ 
sion can be constrained (Piro & Nakar 2013). This is required 
to properly measure the width of the first peak, and thus the 
mass associated with extended material. Spectra taken during 
the first peak are necessary for measuring the energetics and 
composition of the material, so that its radius can be properly 
constrained. Since the range of possible velocities is large, 
relatively low-resolution spectroscopy will be sufficient. 

Further numerical modeling should be conducted to test the 
applicability of these analytic results. This would be useful 
for exploring how efficiently energy is transferred between 
the core and extended material by the traversing shock, for 
which I used Equation (3) here. Such modeling would also 
help explore the impact of the density profile of the extended 
material, whether it be convective, radiative, a wind, or some 
other mass distribution. The recently released open source 
SuperNova Explosion Code (SNEC, Morozova et al. 2015) is 
ideally suited for addressing these questions. 
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thank Giorgos Leloudas, Viktoriya Morozova, Ehud Nakar, 
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